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ABSTRACT 

We present calculations of the spin-down of a neutron star atmosphere due to hydrostatic expansion during a 
Type I X-ray burst. We show that Cumming & Bildsten incorrectly calculated the change in the moment of inertia 
of the atmosphere during a Type I burst, resulting in a factor of two overestimation of the magnitude of the spin- 
down for rigidly rotating atmospheres. We derive the angular momentum conservation law in general relativity, 
both analytically for the case of slow rotation, and numerically for rapidly-rotating stars. We show that contrary 
to the claims of Heyl, but in agreement with Abramowicz and coworkers, general relativity has a small effect on 
the angular momentum conservation law, at the level of 5-10%. We show how to rescale our fiducial results to 
different neutron star masses, rotation rates and equations of state, and present some detailed rotational profiles. 
Comparing our results with recent observations of large frequency shifts in MXB 1658-298 and 4U 1916-053, we 
find that the spin-down expected if the atmosphere rotates rigidly is a factor of two to three less than the observed 
values. If differential rotation is allowed to persist, we find that the upper layers of the atmosphere do spin down 
by an amount comparable to or greater than the observed values. However, there is no compelling reason to 
expect the observed spin frequency to be that of only the outermost layers of the atmosphere. We conclude that 
hydrostatic expansion and angular momentum conservation alone cannot account for the largest frequency shifts 
observed during Type I bursts. 

Subject headings: accretion, accretion disks — nuclear reactions — relativity — stars: neutron — stars: rotation 
— X-rays: bursts 



1. INTRODUCTION 

Type I X-ray bursts are thermonuclear flashes on the sur- 
face of weakly-magnetic (B < 10'" G) accreting neutron stars 
in low mass X-ray binaries (LMXBs) (Lewin, van Paradijs, & 
Taam 1995; Bildsten 1998). Observations with the Rossi X-ray 
Timing Explorer (RXTE) have revealed nearly-coherent oscil- 
lations during Type I X-ray bursts from ten LMXBs, with fre- 
quencies in the range 270-620 Hz (for reviews, see van der Klis 
2000; Strohmayer 2001). These oscillations are interpreted as 
rotational modulation of surface brightness asymmetries, thus 
providing a direct measurement of the neutron star spin fre- 
quency (Strohmayer et al. 1996). The inferred rotation periods 
of a few milliseconds support the idea that the neutron stars 
in LMXBs are the progenitors of the millisecond radio pulsars 
(see Bhattacharya 1995 for a review). 

The large modulation amplitudes, high coherence, and fre- 
quency stability provide convincing evidence that the burst os- 
cillation frequency is related to the neutron star spin frequency 
(Strohmayer 2001 and references therein). However, many puz- 
zles remain. Open questions include (i) the origin of the asym- 
metry on the surface of the neutron star, particularly during the 
fading tail of the burst, when the burning is thought to have 
spread over the whole surface; (ii) why oscillations are seen in 
some bursts but not others, and in some sources but not others 
(e.g., Muno et al. 2001); and (iii) whether the oscillation fre- 
quency is the spin frequency or twice the spin frequency, mo- 
tivated by comparisons of the burst oscillation frequency to the 
difference between the frequencies of the kHz quasi-periodic 
oscillations (QPOs) seen in the persistent emission (van der 



Klis 2000). 

The oscillation frequency increases by a few Hz during the 
burst. To explain this, Strohmayer et al. (1997) proposed that 
the burning shell decouples from the neutron star, and under- 
goes spin changes due to angular momentum conservation as 
it expands and contracts. Cumming & Bildsten (2000, here- 
after CB) studied the rotational evolution of the neutron star 
atmosphere during a burst. Rather than addressing the complex 
problems of how the burning front spreads over the neutron star 
surface during the burst rise, or what causes the asymmetry at 
late times, CB asked what could be learned from detailed mod- 
els of the vertical structure of the atmosphere. They computed 
the vertical hydrostatic expansion during the burst, and the re- 
sulting spin changes. They noted that to see a single, coherent 
frequency requires either the atmosphere be rigidly rotating, or 
the burning region to be vertically thin, so that the differential 
rotation across it is small (see also Miller 2000). In order to 
compare the theoretical results with observations, CB assumed 
that the burning shell rotates rigidly. They found rough agree- 
ment between the calculated spin down and observed values. 

The composition of the burning shell affects the amount of 
expansion. For accretion rates relevant to X-ray bursters, the 
accreted hydrogen (H) burns via the hot CNO cycle as matter 
accumulates on the neutron star surface. The thermally un- 
stable helium (He) burning which leads to the Type I X-ray 
burst may occur either before or after the hydrogen has been 
consumed (see Bildsten 1998 for a review), leading to either 
mixed hydrogen/helium or pure helium ignitions. The presence 
of a substantial amount of hydrogen leads to a factor of two 
greater expansion because of the lower mean molecular weight 
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(/i = 4/3 for pure He; ^ w 0.6 for solar composition). CB noted 
that the largest frequency shifts observed (fractional frequency 
shifts jv w 0.8%) were difficult to achieve with pure He ig- 
nitions, requiring temperatures at the base of the burning shell 
extremely close to the limiting value from radiation pressure at 
the start of the burst. However, they also noted that no burst 
oscillations have been seen during Type I X-ray bursts with 
long 20 s) cooling tails, characteristic of hydrogen burning 
(which is limited by beta decays during the rp process; Bild- 
sten 1998), perhaps indicating a preponderance of He bursts 
amongst those which show oscillations. Thus they suggested 
that understanding the largest frequency shifts observed as be- 
ing due to hydrostatic expansion together with rigid rotation 
might prove problematic. 

Recently, two sources have shown extremely large frequency 
shifts. Wijnands, Strohmayer, & Franco (2001) found an 
increase of « 5 Hz in the 567 Hz burst oscillation from 
MXB 1658-298. In addition to being a large frequency shift 
(Av/v 0.9%), the spin evolution in this burst was unusual, 
with a rapid frequency increase occuring several seconds into 
the burst tail. This frequency shift is a little larger than 
the rigidly-rotating mixed H/He models of CB. Galloway et 
al. (2001) discovered a 270 Hz burst oscillation from the dip- 
ping source 4U 1916-053. During 4 seconds in the burst de- 
cay, the frequency increased by 3.6 Hz, a fractional shift of 
1.3%, by far the largest fractional frequency shift observed so 
far. Comparing with the results of CB, Galloway et al. (2001) 
pointed out that the expansion implied by this frequency shift 
is larger than expected if the atmosphere rotates rigidly, espe- 
cially given that the measured peak flux of the X-ray burst is 
sub-Eddington (w 0.5 Fedd), and that the orbital period of this 
source (w 50 mins; Walter et al. 1982; White & Swank 1982; 
Grindlay et al. 1988) implies a hydrogen-poor companion (Nel- 
son, Rappaport, & Joss 1986) and therefore a small fraction of 
hydrogen present at ignition. 

In this paper, we present new calculations of the hydrostatic 
expansion and resulting spin-down of the neutron star atmo- 
sphere. In §2, we describe an error in CB's calculation of the 
change in the moment of inertia of the burning shell, and show 
that the values of spin-down calculated by CB are a factor of 
two too large. Comparing the new results with observations, 
we show that, even when hydrogen is present at ignition, the 
theoretical values of spin-down assuming rigid rotation are a 
factor of three less than the largest observed frequency shifts. 
We then go on to consider different neutron star masses and 
equations of state. In §3, we derive the angular momentum 
conservation law in general relativity, presenting analytic re- 
sults for slow rotation, and numerical calculations of rapidly- 
rotating stars. We compare our results with recent work by Heyl 
(2000) and Abramowicz, Kluzniak, & Lasota (2001). In §4, we 
show how to include general relativity and rapid rotation in the 
equations describing the structure of the atmosphere. We then 
rescale our fiducial results of §2 to different neutron star masses 
and equations of state, before presenting some detailed models 
and rotational profiles. We again show that the expected spin- 
down due to angular momentum conservation is smaller than 
the observed frequency drifts if the atmosphere rotates rigidly. 
We discuss the implications of our results in §5. 

2. CALCULATION OF MOMENT OF INERTIA AND COMPARISON TO 
OBSERVATIONS 

In this section, we first show that CB made an error of a fac- 
tor of two in their calculation of the change in the moment of 



inertia during a Type I X-ray burst, and then compare the new 
spin down calculations with observations. 

Consider a spherical shell at the surface of the star (The 
arguments of this section are unaffected by including the cen- 
trifugal distortion of the isobaric surfaces.) The mass of the 
shell is 



M-- 



dm ■■ 



47rr p dr, 



and its moment of inertia is 
2 



/ = — / r^ dm = — 
3./ 3 



r p dr, 
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where the factor of 2/3 comes from integrating over angles. 
Following CB, we adopt a plane-parallel approximation for the 
thin atmosphere, and write r = R + z, keeping only first order 
terms in z/R- Changing integration variables from radius r to 
column depth y, where dy = -pdr, equation (^ reduces to CB's 
equation (19) for the moment of inertia. If the angular mo- 
mentum of the atmosphere, lil, is conserved, the spin-down is 
given by comparing its moment of inertia before the burst with 
that during the burst, Ail/D, = -Al/I. 

The error in CB's calculation of the change in moment of in- 
ertia arises because they assume that the total column depth of 
the atmosphere is constant during the burst. We now show that 
this is not the case when the mass of the layer is conserved. We 
write the mass of the shell as 



M = AttR' 



1 + 



2z(y) 
R 



2(z) 

dy = 4TrR^Ay[ 1 + -^ 



(3) 



where Ay is the total column depth, and we define a mass- 
weighted thickness (z) = J z(y)dy/ Ay. If the vertical extent of 
the layer (z) changes, so must its total column depth Ay if the 
total mass M is conserved. Physically, a thin spherical shell has 
a greater surface area if it moves radially outwards, and so must 
have a lower column depth in order to conserve its total mass. 
By taking the column depth at the base during the burst to be 
the same as that before the burst, CB unwittingly increased the 
mass of the layer during the burst, and so overestimated the 
change in its moment of inertia. 

Rather than including the change in the column depth of the 
layer explicitly, we write the mass element in equation (^ as 
dm = AttR^ dy, so that there is a unique correspondence between 
mass and column depth. The moment of inertia is then 



^ttR" 



1 + 



2z(y) 
R 



dy. 



(4) 



where the integration limits are now held fixed during the burst, 
thereby keeping the mass of the layer constant, M = AuR^Ay = 
AnR^iyb-yt)- We have checked the validity of this expression 
by integrating the full spherical equations for this problem. In 
this paper, we adopt a plane-parallel approximation because it 
allows us to straightforwardly incorporate the general relativis- 
tic angular momentum conservation law. 

The spin changes during the burst result from changes in the 
vertical extent of the atmosphere, which are described by the 
second term of equation (^. Equation (19) of CB has a factor 
Az/R rather than 2z/R, so that the spin-down calculated by CB 
is too large by a factor of 2. Specifically, both the spin-down of 
the convection zone shown in Figures 2 and 6 of CB, and the 
spin-down of a rigidly-rotating atmosphere shown in Figure 13 
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of CB, should be reduced by a factor of 2. The spin-down of 
the radiative layers shown in Figures 2, 4 and 6 of CB is, how- 
ever, calculated correctly, since in that case CB assumed r^fl 
was constant for each spherical shell. 

Figure [l| shows the spin-down calculated using equation (Q) 
for the moment of inertia, presuming the whole atmosphere 
rotates rigidly during the burst. This is an updated version 
of Figure 13 of CB, and detailed discussion of the theoreti- 
cal models may be found in that paper The calculated spin 
down is exactly a factor of two less than found by CB. We take 
g= \.9 X 10'"* cm s"^, the Newtonian gravity for a 1.4 Mq, 10 
km star. We consider two pre-burst models, mixed H/He ig- 
nition at m = 0.1 iriEdd (solid lines and squares), and pure He 
ignition at m = 0.015 iriEdd (dotted lines and open squares). 

Figure [l| shows results for both convective and radiative mod- 
els. For the radiative models, we show the spin down as a func- 
tion of flux, since the radiative solution is uniquely determined 
by the flux. The convective models depend both on the flux and 
the extent of the convection zone; we plot the spin down result- 
ing from the convective models given in Table 3 of CB. The 
convective models have a composition profile the same as that 
prior to the burst. The radiative models have either a pre-burst 
composition (upper solid and dotted curves) or a composition of 
heavy elements (lower solid and dotted curves), either ^^Ni for 
pure He burning, or ^^Kr for mixed H/He burning (we choose 
^^Kr as representative of the products of rp process H burning; 
see Schatz et al. 1998; Koike et al. 1999; Schatz et al. 2001). 
As described in CB, the radiative and convective models with a 
pre-burst composition are appropriate for the burst rise, whereas 
the radiative atmospheres composed of heavy elements are ap- 
propriate for the cooling tail of bursts. Convective models spin 
down more than radiative models because the steeper tempera- 
ture gradient in the convection zone leads to a larger base tem- 
perature for a given flux. 

On the left of Figure 1, we plot the observed frequency drifts 
as horizontal bars. A summary of these observations may be 
found in CB, Table 1; in addition, we include the frequency 
shifts recently reported for 4U 1916-053 (Gafloway et al. 2001) 
and MXB 1658-298 (Wijnands et al. 2001). For each source, 
we plot the largest observed frequency shift. We plot as upper 
limits those cases for which a frequency was only seen in the 
burst tail. The largest observed frequency shift is larger than the 
theoretical values by at least a factor of 3. 

3. ANGULAR MOMENTUM CONSERVATION IN GENERAL 
RELATIVITY 

In this section, we derive the angular momentum conserva- 
tion law in general relativity. Our approach is to calculate the 
conserved angular momentum of a particle moving on a circular 
path (in the Appendix, we show that the same result is obtained 
starting with the conservation equations for a fluid). We show 
analytic results for slowly rotating stars, and present numeri- 
cal calculations of rapidly rotating stars. We compare our re- 
sults with recent calculations by Heyl (2000) and Abramowicz 
et al. (2001). 

In Newtonian physics, the specific angular momentum of a 
particle is given by r^fl, where r is the distance from the ro- 
tation axis. If the particle's motion conserves angular momen- 
tum, the change in its spin frequency with radius is given by 
dlnil/dlnr = -2. In this section, we derive the equivalent ex- 
pression in general relativity. We start by discussing the angular 
momentum of a particle moving in the Schwarzschild metric. 
We then derive the conserved angular momentum in the slow 



rotation approximation, i.e. neglecting terms of order (n/Qn)^, 
where fi^ = GM/R^, allowing us to write an analytic expression 
for the angular momentum conservation law. However, since 
(n/riK)^ may be quite large (for example, {H/Hk)^ = 0.26 for 
a star with 7? = 15 km and M = 1.4M0 spinning at 600 Hz), 
we conclude this section with detailed numerical calculations 
of rapidly rotating models. 

Consider first a spherical star, for which the metric is 



ds^ = -e''dt^ + e^dr^ + r^ \ 



(5) 



where outside the star the metric functions A and v are given by 



1- 



2M 



(6) 



where M is the total gravitational mass. Note that in this paper 
we write the ratio GM /rc^ asM/r for clarity, but otherwise ex- 
plicitly write in the constants G and c. A particle moving in a 
circular path has four-velocity (e.g., Hartle 1970) 



(7) 



where 0(r) is the angular velocity measured by an observer at 
infinity, u' is the redshift factor between the particle and the 
observer at infinity, and f" and 0" are Killing vectors of the 
spacetime. The redshift factor is 



-1//2 . 



^_2My'/2 



(8) 



where we choose the standard normalization gapu^u^ = -1. 
Associated with the Killing vectors f " and are the conserved 
quantities E = -u^ta and L = u"ipa, the specific energy and an- 
gular momentum. Using equations (^ and we obtain (see, 
e.g. Schutz 1990, §7.4) 



and 



2M 



L = r^sin^enl 1 



1/2 



2M 



-1/2 



(9) 



(10) 



Differentiating equation ( pX| ) and evaluating it at the surface of 
the star, we find 



dlnft _ f 5M 
dlnr ~~ V ~2^ 



2M 



(11) 



when a particle moves such that angular momentum is con- 
served. In the limit M <^R, this formula reduces to the Newto- 
nian result. 

It is simple to extend this analysis to the case of slow rotation, 
for which the metric is (Hartle 1967) 



£// = -e^df" + e^dr-" + AdO"- + sin^ Odcf)^) - 2ujr^ sin^ 9 d(t)dt, 

(12) 

where A and ly have the same solutions as the spherical case 
outside the star, and where uj = 2J /P outside the star (7 is the 
total angular momentum of the star). Again solving for the con- 
served angular momentum, we find 



L = ?5m^e {n-Lj){ 1-— 



-1/2 



(13) 
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giving 



dlnr 



5M I ( , M 

1 + —r 1 

2R R 



.-^1 , 04) 



at the surface of the star, where we presume the particle is ini- 
tially rotating with the spin of the star, and write the moment of 
inertia I = J /Vl. 

For the rapidly-rotating case, we consider the general metric 
for an axisymmetric stationary star, 

ds^ =-e^^Pdf + e''-Pr'^?,\n^e{d(j}-u}dtf + e^°' {dr^ + fdO^) , 

(15) 

where the metric potentials p^^^a and lo depend only on 
the coordinates r and 6. The coordinate F is related to the 
Schwarzschild coordinate r by r = rexp(|(7 - p)), so that 
InrsmO is the circumference of a circle with constant r and 
9. For further details about the interpretation of these coordi- 
nates and potentials, see Friedman, Ipser & Parker (1986) and 
Morsink & Stella (1999). Given an equation of state, a mass and 
a rotation rate, the metric (|l5]) can be solved numerically. We 
use a code written by N. Stergioulas'*, which assumes rigid rota- 
tion and is based on methods developed by Komatsu, Eriguchi 
& Hachisu (1989) and Cook, Shapiro & Teukolsky (1994). 

We present the properties of some rapidly-rotating models in 
Table 1. We consider two different equations of state (EOS), 
EOS L which is stiff and EOS APR which is softer. EOS L 
(Pandharipande & Smith 1975) is one of the stiffest EOS in 
the Arnett & Bowers (1977) catalogue. EOS APR is the model 
A18h-(5vh-UIX* computed by Akmal, Pandharipande & Raven- 
hall (1998) which uses modern nucleon scattering data and first 
order special relativistic corrections. For each EOS, we con- 
sider two different masses, M = 1.4 Mq and 2.0 M0, and two 
spin frequencies, v = 300 Hz and 600 Hz. The last column in 
Table 1 gives the break-up spin frequency in each case. 

We work in the equatorial plane, and do not consider the vari- 
ations of the metric potentials with latitude. The specific angu- 
lar momentum is 

vr 

L=^=, (16) 



(see also Morsink & Stella 1999) where the three-velocity of a 
corotating particle as measured by a zero angular momentum 
observer is given by 



-(7+p)/2 



[n{r)-Lu(r)]. 



(17) 



Again holding L constant as we vary the Schwarzschild coordi- 
nate radius r, we find 



dlnn 
dlnr 



1-y- 4 (7r+Pr) 



(18) 



where R is the equatorial radius of the star, subscripts denote 
partial differentiation and all quantities are evaluated at the 
equator. In the slow rotation limit, equations ( p^ and ( p^ re- 
duce to equations ( pj[ ) and (]I4|). 

Heyl (2000) recently calculated In fi/c/ In r, with a rather 
different result (compare our eq. | ]l4| ] with his eq. [10]). How- 
ever, Abramowicz et al. (2001) point out a sign error in Heyl's 
calculation, and, more importantly, that Heyl assumes that the 
quantity L/E is conserved, as, for example, a particle in orbit 
(Abramowicz & Prasanna 1990). However, the energy E of 



a fluid element in the atmosphere is not conserved during the 
burst. The correct conserved quantity, which we have consid- 
ered in this section, is the angular momentum per particle L. 
In the slow rotation approximation, our result (eq. [|l4||) agrees 
with equation (9) of Abramowicz et al. (2001). 
We write dlnfl/dlnr as -2/3, where 



2 dlnr 



(19) 



is unity in the Newtonian limit. We give the value of /3 in Table 
[1] Figure^ shows /3 as a function of neutron star mass for EOS 
APR and EOS L, and for ly = 300 and 600 Hz. We see that in- 
cluding the general relativistic angular momentum conservation 
law changes dlnft/dlnrhy 5-10%. 

4. CALCULATION OF EXPANSION AND SPIN DOWN 

In this section, we present calculations of the expansion and 
spin-down of the atmosphere, including the effects of general 
relativity and rapid rotation. In the spirit of CB, we ignore lat- 
itudinal variations, and work in the equatorial plane. First, in 
§4.1, we discuss the effects of general relativity on the equa- 
tions describing the hydrostatic and thermal structure of the at- 
mosphere, and calculate the reduction in gravity at the equator 
due to rapid rotation. In §4.2, we show how to include the gen- 
eral relativistic angular momentum conservation law that we 
obtained in §3. In §4.3, we rescale our fiducial results of §2 
for the spin-down of a rigidly-rotating atmosphere to different 
masses, equations of state, and spin frequencies. In §4.4, we re- 
lax the assumption of rigid-rotation, and present the rotational 
profiles of the atmosphere for several different cases. 

4. 1 . Hydrostatic and Thermal Structure of the Atmosphere 

We now write down the differential equations describing the 
hydrostatic and thermal structure of the atmosphere including 
general relativity, and compare them to the Newtonian equa- 
tions. We start by considering a non-rotating star, and then dis- 
cuss the additional effects of rotation. 

For a spherical star, Thorne (1977) (see also Thorne 1967) 
gives the equation for mass conservation as 



dMr 



= 4Trr^pV{r), 



(20) 



where M,- is the rest mass (number of baryons multiplied by 
baryon rest mass) within coordinate radius r, p is the rest mass 
density, and V is the "volume redshift factor". For slow rota- 
tion, V(r) = (1 —2M{r)/r)~^l^, where M(r) is the gravitational 
mass interior to r. We adopt a plane parallel approximation 
in the thin envelope, and write r = R + z/V, where z ^ is 
the proper length. Here, V is the redshift factor at the surface 
of the star V = V(R). We define the rest mass column depth 
dy = -dMr/4-TTR^. The mass conservation equation becomes 



dz 
d'y 



The equation for hydrostatic balance is (Thorne 1977) 
dP GM(r)^ 



dMr 



47rr4 



■V, 



(21) 



(22) 



Code publicly available at http : / / www . gravity . phys . uwm . edu/rns / 
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where P is the locally measured pressure, and we neglect contri- 
butions to the gravitational mass from the fluid energy density. 
Using the definition of column depth, we rewrite equation (g^) 
for the thin layer as 

<^'> 

dy 



where 



GM 



V. 



(24) 



Rewriting the entropy and flux equations (Thorne 1977, eqs. 
[lid] and [llh]) in a similar way, we find 



dF 

dy 



and 



dT 
dy 



3k 
AacT^ 



(25) 



(26) 



where F is the heat flux, e is the nuclear energy release rate 
per unit rest mass, n is the opacity, and T is the locally mea- 
sured temperature. Equations (pl|), (^3|), (^^, and (|2^, which 
describe the structure of the atmosphere, are identical to the 
Newtonian equations (compare CB §2). 

These results also hold in the slow rotation approximation, 
for which the star is spherical to first order in (Vl/VLKf'. The in- 
clusion of rapid rotation has three effects. The first is to change 
the volume correction factor V. In addition, the mass conserva- 
tion equation must be covariant so that all observers agree on 
the rest mass of the star This requirement gives an additional 
Lorentz factor in equation ( |20| ) (e.g., Friedman et al. 1986). 
However, we find that for the models we consider in this pa- 
per, both of these effects are < 1%, and so we neglect them. 

Much more important is the correction to gravity due to the 
rapid rotation. This correction factor is the relativistic general- 
isation of the Newtonian reduction in gravitational acceleration 
due to the centrifugal force. To calculate this, we write the 
four-velocity as m" = ??.(?" + ili/)"), and use the normalization 
m^Mq = -1 to solve for the redshift factor 



1 



(27) 



where v is given by equation (|l7|). The equation of hydrostatic 
balance for a rigidly rotating star is then (Friedman et al 1986) 



IdP d ^ ^ 
-— =--log7^, 
p dr dr 



(28) 



where we assume the pressure is much smaller than the rest 
mass density. To reduce this equation to the form of equation 
(^3|), we write 

.= ^VA (29) 

where the correction factor A accounts for the reduction in g at 
the equator due to rapid rotation. 



A = 



GMV^ dr 



log7^, 



(30) 



tending to unity in the limit of zero rotation. This is the most 
important effect coming from rotation, giving a correction of up 
to 25% for a rotation frequency of 600 Hz and a stiff equation 
of state. 



In summary, the general relativistic equations look the same 
as the Newtonian equations if we use the correct value of sur- 
face gravity g (eq. p^]), and if we identify the column depth y 
with the rest mass column depth, and the thickness z with the 
proper thickness of the layer. Values of g, V, and A are given 
in Table |l] for different models. 

4.2. Calculation of Spin Down 

To incorporate the general relativistic angular momentum 
conservation law into our calculations of spin down, we write 
the angular momentum per unit rest mass at the equator in the 
form 

u^ = mr), (31) 

where /(r) is a function of radius, e.g. /(r) = in the Newto- 
nian limit. If the angular momentum is constant, we find that 
the change of with radius is determined by the gradient of 
/ with radius: [3 = (-l/2)(dlnn/dlnr) = {l/2)(dlnf/dlnr). 
Now consider a particle which moves from radius RtoR + Ar, 
where Ar ^ R. Conservation of angular momentum gives 
u^{R + Ar) = u^{R). Expanding in a Taylor series around r = R, 
we find 

In the Newtonian limit, for which (3= 1, we recover the familiar 
result Ari /VL = -2/S.r/R. To calculate the moment of inertia, we 
must modify equation (Q). A gain writing the angular momen- 
tum in the form of equation (IT)), and expanding f{r) around 
r = R,we find 



/ = 



^TTR^m 



1+ 



2pz(y) 



VR 



dy, 



(33) 



which reduces to equation (m in the Newtonian limit. 



4.3. Scaling our Fiducial Results to Different Masses and 
Equations of State 

The calculations presented in §2 were for a surface grav- 
ity g = 1.9 X 10'^ cm s"^ (the Newtonian gravity for a 1.4M0, 
R=10 km neutron star). We now rescale these fiducial results 
to different masses and equations of state, including the effects 
of rapid rotation and general relativity. 

To calculate the thickness of the atmosphere in §2, we in- 
tegrated hydrostatic balance in the form dP/dz = -pGM /R^, 
giving the thickness /S.z/R oc 1 /gR oc R/M (see also §2 of CB). 
In addition, we must include a factor V to obtain the correct 
general relativistic value for g, and a factor A to include the ef- 
fects of rotation on gravity. We thus find the fractional change 
in thickness of the atmosphere is 



R \M/R J 



Az 
R 



(34) 



where {Az/R){) is the fiducial value from §2. The fractional 
change in spin frequency is f3(2Ar/R), where Ar is the coor- 
dinate thickness. Converting to proper thickness Az involves 
another factor of V, giving 



n [m/r 



/An\ 



An\ 



(35) 
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where we define a as the scaUng factor from the fiducial value, 
and we took M/R = 0.21 in §2. Note that the largest contribu- 
tions to a come from rescaling gravity and including the cen- 
trifugal force, rather than from general relativistic effects. 

Figure H shows a as a function of neutron star mass, for EOS 
APR and EOS L, and in each case for ly = 300 and 600 Hz. We 
also give a in Table |l} Figure |] shows that for massive neutron 
stars with M w 2Mq, the frequency shifts calculated in §2 and 
shown in Figure 1 should be multiplied by 0.65-0.85 for EOS 
L, and 0.3-0.4 for EOS APR. For a IAMq star, these factors 
are 1.1-1.7 for EOS L, and 0.7-0.8 for EOS APR. The largest 
values of a are obtained for a low neutron star mass, stiff equa- 
tion of state, and rapid rotation. However, even a = 1.7 for 
M = 1.4 Mq, EOS L, and u = 600 Hz is not lai-ge enough to 
give agreement with the largest observed frequency shifts. 

4.4. Rotational Profiles 

In Figure ^ we present some rotational profiles for some spe- 
cific models. We take the neutron star mass to be 1.4 Mq, the 
global accretion rate to be 0.1 Mndd, and again work in the 
equatorial plane. We assume that the atmosphere is rigidly- 
rotating immediately prior to the burst, and that during the 
burst, the atmosphere is radiative and carries a flux equal to 
the solar-composition Eddington flux at the photosphere (/^Edd = 
8.8 X 10^4 gi4 erg cm'^ s"'). For ly = 300 and i/ = 600 Hz, and 
for EOS L and EOS APR, we plot the rotational profiles assum- 
ing either (i) complete rotational coupling, resulting in rigid ro- 
tation across the layer, or (ii) no angular momentum transport, 
giving a differentially-rotating layer. 

If we allow differential rotation to persist, we see that the up- 
per layers of the atmosphere spin down by an amount compara- 
ble to or greater than the observed spin changes during bursts. 
However, it is not clear why the spin frequency observed should 
be that of only the outermost layers. Indeed, CB (§3.4) argued 
that substantial differential rotation would wash out any sig- 
nal from the deeper cooling layers because of the finite time to 
transport heat vertically. 

5. SUMMARY AND DISCUSSION 

We have presented new calculations of the hydrostatic ex- 
pansion and spin-down of a neutron star atmosphere during a 
Type I X-ray burst. Our main conclusion is that hydrostatic ex- 
pansion is not enough to explain the observed frequency drifts 
during Type I X-ray bursts if the burning atmosphere rotates 
rigidly. 

In §2, we showed that Gumming & Bildsten (2000) (CB) 
overestimated the change in the moment of inertia of the atmo- 
sphere, obtaining values of spin-down that were a factor of two 
too large. Figure 1 compares the new calculations of spin-down 
with observations, including recent measurements of large fre- 
quency drifts during bursts (Galloway et al. 2000; Wijnands et 
al. 2001). We find that the largest observed frequency shift is a 
factor of 3 or more greater than the theoretical values. 

In §3, we derived the angular momentum conservation law 
in general relativity. We calculated the variation of spin fre- 
quency with radial distance, dlnil/dlnr, for a particle mov- 
ing with constant angular momentum. In the slow rotation 
approximation (fl^ ^ GM/R^), we obtained the analytic re- 
sult given by equation (p^, which agrees with recent work by 
Abramowicz et al. (2001). For rapidly-rotating stars, we cal- 
culated (ilnil/iilnr by numerically solving for the structure of 
the neutron star The correction to the Newtonian angular mo- 
mentum conservation law, (3 = (-1 /2)(dlnil/dlnr), is shown in 



Figure 2 for different neutron star masses, equations of state and 
spin frequencies. Contrary to the results of Heyl (2000), which 
were also shown to be incorrect by Abramowicz et al. (2001), 
we find that the general relativistic correction is small, about 
5-10%. 

In §4, we calculated the atmospheric expansion and spin- 
down, including the effects of rapid rotation and general rel- 
ativity. Working in the equatorial plane (in the spirit of CB's 
calculation, we neglect latitudinal variations in this paper), we 
calculated the scaling factor a required to rescale our fiducial 
results of §2 to different neutron star masses, equations of state, 
and spin frequencies (Figure 3). In addition, we presented the 
rotational profiles for some particular models (Figure 4). We 
find that the largest spin-down is for rapid rotation, and low 
mass stars with a stiff equation of state. For example, a 1 .4 Mq 
star spinning at 600 Hz with the stiff equation of state EOS L 
has a =1.7. However, this is not a large enough factor to bring 
observations and theory into agreement in Figure 1 . 

The rotational profiles given in Figure 4 show that in princi- 
ple frequency shifts as large as those observed can be obtained 
by hydrostatic expansion, if we consider only the outermost lay- 
ers of the atmosphere, and allow differential rotation. However, 
it is not at all obvious why the observed frequency would be 
that of the outermost shells of the atmosphere, which contain 
a small amount of the mass, particularly since the energy re- 
lease is in the deeper layers. Indeed, CB argued (see their §3.4) 
that substantial differential rotation in the burning layers would 
wash out the signal, because of the finite time needed to trans- 
port heat vertically in the atmosphere. 

Another possible objection to the angular momentum con- 
servation picture was given by CB, who pointed out that if the 
burning layers are threaded by a large scale poloidal magnetic 
field, this field will be wound up by the differential rotation 
during the burst. The wound up toroidal field acts back on 
the shear, halting and reversing its direction in the time for an 
Alfven wave to cross the atmosphere (see, for example. Spruit 
1999). For a magnetic field typical of a millisecond radio pul- 
sar, B ~ 10^ G, this timescale is only ~ 0.01 s. One possibility 
is that the surface field is much weaker, B <J 10^ G, allowing 
shearing to persist. Observations of burst oscillations in the 
accreting millisecond X-ray pulsar SAX J1808. 4-3658, which 
perhaps has a - 10**-10'' G field (Psaltis & Chakrabarty 1999), 
would give an interesting test of this picture. In't Zand et al. 
(2000) report a marginal detection with BeppoSAX of oscil- 
lations at 400 ± 2 Hz during a Type I burst, but are unable to 
resolve any frequency drift. Unfortunately, as yet no burst os- 
cillations have been detected from this source with RXTE. 

The results we have obtained in this paper suggest that we 
may have to look elsewhere for an explanation of the frequency 
drifts. Spitkovsky, Levin, & Ushomirsky (2001), in a detailed 
investigation of hydrodynamic flows during Type I bursts, pro- 
pose that the frequency drifts may be explained by a combina- 
tion of radial expansion and variations in the velocity of zonal 
flows on the neutron star surface. Another possibility men- 
tioned by CB is that the burst oscillation is due to a non-radial 
oscillation in the neutron star surface layers. Such possibilities 
remain to be investigated in detail. 
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Table 1 

Neutron Star Models 



M(Mq) 


/?(km) 


M/R 


V" g (1014 cm 






a" 






EOS APR V = 300 Hz 


1.4 
2.0 


11.5 
11.1 


0.18 
0.27 


1.25 1.72 
1.46 3.10 


0.98 
0.99 


0.94 
0.90 


0.72 
0.34 


1070 
1330 


EOS APR I' = 600 Hz 


1.4 
2.0 


11.9 
11.3 


0.17 
0.26 


1.24 1.48 
1.45 2.83 


0.91 
0.95 


0.94 
0.90 


0.81 
0.37 


1070 
1330 


EOS L 1/ = 300 Hz 


1.4 
2.0 


15.2 
15.3 


0.14 
0.19 


1.17 0.90 
1.28 1.39 


0.95 
0.97 


0.96 
0.94 


1.13 
0.65 


730 
850 


EOS L = 600 Hz 


1.4 
2.0 


16.7 
16.2 


0.12 
0.18 


1.15 0.55 
1.25 1.05 


0.72 
0.83 


0.95 
0.93 


1.70 
0.82 


730 
850 



''The redshift factor V = (1 -2M/R) (we neglect the ^1% correction due to rotation). 

''The correction to gravity due to rapid rotation, g = GMVA/R^. 

'^The general relativistic correction to the angular momentum conservation law, (3 = 
i-l/2){d\nCt/d\nr). 

"^The multiplicative scaling factor that should be applied to the results of §2 and Fig. 1, 
to allow for the different gravity, and the effects of general relativity and rapid rotation (see 
text). 

^The break-up spin frequency. 
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n I I rn^ I I I 1 1 1 rn^ 
4U 1916-053 



■ H/He 
□ Pure He 




Fig. 1. — The spin evolution of the atmosphere for convective (squares) and radiative models (lines), assuming rigid rotation is maintained throughout the 
atmosphere. We plot the fractional spin frequency change of the burning shell, AQ/U, such that AQ/Q > indicates spin down. We show results for both mixed 
H/He ignitions (solid lines, solid squares) and pure He ignitions (dotted lines, open squares). For the convective models, the composition is the same as before 
ignition. For the radiative models, the upper curve is for a pre-burst composition; the lower curve is for a composition of ^*Ni (pure He) or ^''Kr (mixed H/He). We 
indicate the observed frequency shifts by horizontal bars. These observations are summarized in Table 1 of CB; in addition, we include the frequency shifts recently 
reported for 4U 1916-053 (Galloway et al. 2001) and MXB 1658-298 (Wijnands et al. 2001). For those bursts in which the oscillation frequency was seen only in 
the tail, we plot the AQ/Q value as a lower limit. This Figure is an update of Figure 1 3 of CB . 
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— EOS APR, 300 Hz - 

— EOS APR, 600 Hz _ 
--EOS L, 300 Hz 

— - EOS L, 600 Hz 




0.85 



1.2 



1.4 



I I I I I 
1.6 1.8 
Mass (Mq) 



J L 



J L 



2.2 



Fig. 2. — The general relativistic correction to the Newtonian angular momentum conservation law. We plot /3 = (-l/2)(dlnQ/dlnr) as a function of neutron star 
mass for EOS L (dashed curves) and EOS APR (solid curves), and for u = 300 Hz (light curves) and u = 600 Hz (heavy curves). 
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Fig. 3. — The scaling factor a which rescales the fiducial values of AQ/Q presented in §2 and Figure 1 to different neutron star masses, equations of state, 
and spin frequencies. The factor a includes the scaUng with surface gravity, the correction to gravity at the equator due to rapid rotation, and general relativistic 
corrections to gravity and the angular momentum conservation law. We show results for EOS L (dashed curves) and EOS APR (solid curves), and for v = 300 Hz 
(light curves) and v = 600 Hz (heavy curves). 
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Fig. 4. — Detailed rotational profiles for 1.4 Mq neutron star.s undergoing mixed H/He ignition at a global accretion rate of 0.1 MEdd- We as.sume that the 
atmosphere is rigidly-rotating immediately prior to the burst, and that during the burst, the atmosphere is radiative and carries a flux equal to the solar-composition 
Eddington flux at the photosphere. For v = 300 Hz (left panel) and u = 600 Hz (right panel), and for EOS L (solid lines) and EOS APR (dotted lines), we plot 
the rotational profiles assuming either (i) complete rotational coupling, resulting in rigid rotation across the layer, or (ii) no angular momentum transport, giving a 
differentially-rotating layer. 
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APPENDIX 



ANGULAR MOMENTUM CONSERVATION FOR A FLUID 

In this Appendix, we start with the conservation equations for a fluid, which include contributions from the fluid's enthalpy and 
from heat flow, and show that the appropriate conserved quantity is the angular momentum per baryon as calculated in §3. 
We begin with the stress tensor for a fluid 



(see, for example, Mihalas & Mihalas 1984). The enthalpy h is 

e+P 



h = 



P 



(Al) 



(A2) 



where the energy density 

e = pc^ + p\l-pB (A3) 

is the sum of rest mass energy, internal energy, and nuclear binding energy. We write the heat flux as F". By setting ^^Vq-T"'^ = 0, 
and using baryon number conservation Vq(/9m") = 0, we obtain the angular momentum conservation law 



pM^Vc, [hu^+^^+Va (f"M0) =0. 



(A4) 



Assuming the azimuthal heat flux F"^ may be neglected when compared to the vertical heat flux F'' = F, we rewrite equation ( A4 ) as 

(A5) 



4" (hu^) =---^ (Fua,) . 
dt ^ ^' pdz^ ' 



We now simplify this equation using the first law of thermodynamics, which is 



dt ■■ 



P 



dp + pTds 



where s is the entropy per baryon mass, or, using equation ( A2), 



dP 

dh = Tds H . 

P 



Substituting for dh/dt in equation ( A5), we find 

du^ 
dt 



ds 1 dF 

T — + 

dt p dz 



dP F du^ 
ph dt ph dz 



(A6) 



(A7) 



(A8) 



The first term in equation ( AS ) vanishes, since Tds /dt = -{ \/ p)dF /dz is the familiar entropy equation. The nuclear energy production 
rate is included in the Tdsjdt term, since the nuclear binding energy is included in the enthalpy: equations (A3) and ( A6) give 



^ds _dli P dp dB 
dt dt p^ dt dt ' 



(A9) 



where dB/dt is the nuclear energy generation rate. Integrating the second term in equation ( \8 ) in time, we find it gives a contribution 
to Am0/m0 of order AP/pc^ (H/R)(P/pc^) « (gR/c^)(H/Rf w 10"^ where we take the pressure scale height « 10 m. The 
third term in equation ( |A^ ) represents radiative viscosity (see Mihalas & Mihalas 1984). Again integrating in time, we find that this 
term gives a contribution sa {F / pc^){t / H){/S.u^/ U(j,) w lO~'*(f/sec)(AM0/M0). Both the second and third terms are of order (H/R)^ 
or smaller, and may be neglected for our purposes. 

In summary, consideration of the conservation law for a fluid gives 



du^ 

~dr 



= 0+0 



(AlO) 



showing that, at the level of accuracy we require, the angular momentum per baryon is conserved. 



